Spinon deconfinement in doped frustrated quantum antiferromagnets 
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The confinement of a spinon liberated by doping two-dimensional frustrated quantum antiferro- 
magnets with a non-magnetic impurity or a mobile hole is investigated. For a static vacancy, an 
intermediate behavior between complete deconfinement (kagome) and strong confinement (checker- 
board) is identified in the J1 — J2 — J3 model on the square lattice, with the emergence of two length 
scales, a spinon confinement length larger than the magnetic correlation length. For mobile holes, 
this translates into an extended spinon-holon boundstate allowing one to bridge momentum (ARPES 
spectral function) and real space (STM) experimental observations. These features provide clear 
evidence for a nearby "deconfined critical point" in a doped microscopic model. 

PACS numbers: 75.10.-b, 75.10.Jm, 75.40.Mg 



The search for exotic spin liquids (SL) has been enor- 
mously amplified after the discovery of the high criti- 
cal temperature (high-Tc) cuprate superconductors. In- 
deed, Anderson suggested that the Resonating Valence 
Bond state is the relevant insulating parent state that 
would become immediately superconducting under hole 
doping ^J. Such a state is characterized by short range 
magnetic correlations and no continuous (spin) or dis- 
crete (lattice) broken symmetry. Another major charac- 
teristic of this SL phase is the deconfinement 2] of the 
S=l/2 excitations (spinons) in contrast to ordered mag- 
nets which have 5 = 1 spin waves. Upon doping, some 
scenarii predict a 2D Luttinger liquid 3], i.e. a state 
which exhibits spin-charge separation, a feature generic 
of one-dimensional correlated conductors. 

Magnetic frustration is believed to be the major tool to 
drive a two-dimensional (2D) quantum antiferromagnet 
(AF) into exotic quantum disordered phases. The Va- 
lence Bond Solid (VBS), an alternative class of quantum 
disordered phases breaking lattice symmetry, seems to be 
a strong candidate in some frustrated quantum magnets 
as suggested by robust field theoretical arguments |4|, 
early numerical computations of frustrated quantum AF 
on the square lattice with diagonal bonds @, |y, an d m 
the 2D checkerboard lattice |S| (with diagonal bonds only 
on half of the plaquettes). In contrast, the 2D kag ome 
lattice 9] shows no sign of ordering of any kind |ld Tll| . 
The "deconfined critical point" (DCP), a new class of 
quantum criticality, was proposed recently e.g. to de- 
scribe a continuous AF to VBS transition [l2|. 

Investigation of hole doping in frustrated magnets 
has revealed stricking differences between VBS and SL 
phases although they both exhibit a finite spin-spin cor- 
relation length £af- Viewing these phases as fluctuat- 
ing singlet backgrounds, removing an electron at a given 
site (e.g. by chemical substitution with an inert atom) 
or, as in Angular Resolved Photoemission Spectroscopy 
(ARPES) experiments, in a Bloch state of given momen- 
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FIG. 1: (Color on-line) Schematic picture of a vacancy (or 
doped hole) in a frustrated magnet. The segments stand for 
singlet bonds and the arrow represents the spinon liberated in 
the process, (a) Holon-spinon BS in a columnar VBS bound 
by a "string" potential (dotted line), (b) Deconfined holon 
and spinon in an (hypothetical) SL host. 



turn naturally breaks a spin dimer and liberates a spinon, 
i.e. a S=l/2 polarisation in the vicinity of the empty site 
(holon). The single hole spectral function shows a sharp 
peak (resp. a broad feature) characteristic of a holon- 
spinon boundstate (resp. holon-spinon scattering states) 
in the checkerboard VBS phase (resp. kagome SL phase). 
The simple physical pictures behind these two typical be- 
haviors are depicted in Fig. ^a) and (b) for a confining 
columnar dimer phase and a SL phase respectively. The 
new length scale £ con f (average distance between vacancy 
and spinon) which emerges naturally in the VBS phase 
is to be identified with the correlation length over which 
dimer (or VBS) order sets in. Interestingly, it has been 
predicted that, in the vicinity of the DCP, confinement 
occurs on a much larger length scale £ con f which diverges 
as a power law of £af • 

So far the DCP scenario is only supported by field the- 
oretic arguments. It is therefore of crutial importance to 
investigate its relevance in the framework of microscopic 
models. In this Letter, we address the issues depicted in 
Fig. [l] by considering a single hole introduced in the 2D 
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spin- 1/2 AF J1-J2-J3 Heisenberg model on the square 
lattice at zero temperature defined by 

H = J2JijSi-$j (1) 

(ij) 

where the exchange parameters are limited to first 
(Ji), second (J 2 ) and third (J 3 ) N.N. AF couplings. The 
classical phase diagram of this model [30 is very rich 
(see Fig. |2| showing four ordered states - Neel, collinear 
(q = (tt, 0)) and two helicoidal - separated by contin- 
uous or discontinuous boundaries. The subtle interplay 
between quantum fluctuations and frustration (J2 and 
J3 terms) is expected to destabilize the classical phases 
and lead to a quantum disordered singlet ground state, 
possibly of VBS type. We show that one of the major 
prediction of the DCP scenario, namely the emergence 
of a hierarchy of length scales is indeed observed for in- 
termediate frustration in correlation with the possibility 
of a direct Neel- VBS continuous transition. This finding 
is contrasted to two other extreme behaviors -complete 
deconfinement and strong confinement- observed in the 
kagome and checkerboard lattices respectively. 

Let us first briefly review some results in the literature 
supporting the existence of a cristaline quantum disor- 
dered phase in model {1} (leading to spinon confinement). 
In the parameter range where frustration is largest, many 
approaches, including spin- wave theory ^(|, exact diag- 
onalizations || , series-expansion and large- TV expan- 
sions 2], have firmly established for J3 = the relative 
stability of a quantum disordered singlet ground state: 
a columnar valence bond solid with both translational 
and rotational broken symmetries Q or a plaquette state 
with no broken rotational symmetry have been pro- 
posed. For the pure J\ — J3 model, a non-classical phase 
also appears between the Neel (jr. tt) and the spiral (g, q) 
phases : a VBS columnar state 1 181 or a succession of a 
VBS and Z2 spin-liquid phases have been proposed. 
Lastly, when J 2 and J3 are both non-zero, in the range 
(J3 + J2)/ J\ ~ 0.4 — 0.6, the finite size scaling anal ysis 
of the dimer susceptibility computed up to 50 sites |2£| 
shows a non- vanishing signal again strongly suggesting a 
VBS order Q. 

A static vacancy at a given site O of the lattice is 
a simple setup relevant to test the DCP ideas and to 
experiments. In pratice, the vacancy is simply modelled 
by setting to zero all the couplings involving site O 
and the computations are performed by Lanczos ED of 
a cluster of 32 sites (i.e. V^32 x a/32) which respects all 
point group symmetries of the infinite lattice. Such an 
impurity acts, theoretically, as a local probe of the host. 
It can be viewed alternatively as a localized holon (S = 
and charge Q — e) so that the form of the surrounding 
spin density is expected to provide valuable insights on 
the spin-charge confinement /deconfinement mechanism. 

The single impurity Green function G(uj) = 
(^bareK^ — ^) _1 |^bare) is computed by (i) constructing 
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FIG. 2: (Color on-line) Classical phase diagram for the J\- 
J2-J3 model. Second order (discontinuous) transitions are 
indicated by dashed (solid) lines (see e.g. Ref. H3)- The 
shaded (blue online) region shows the approximate location 
of the minimum of the spectral weight Z in the quantum 
version. The region with a weight between 0.79 and 0.84 on 
the 32-site cluster is delimited by dashed lines and red dots. 



the (normalized) "bare" initial state |^bare) = 2co,cr|^o) 
from the host GS \^o) by removing an electron of spin 
a and (ii) using a standard Lanczos continued-fraction 
technique. Most of the cj-integrated spectral weight (nor- 
malized to 1) of ImG(w) is in fact contained in the lowest 
energy pole of weight Z = |(^ gs |^bare)| 2 where |^ gs ) is 
the (normalized) GS of the system with one vacancy at 
site O. Results shown in Fig. |H{a,b) show however that 
Z is significantly suppressed in the region where a quan- 
tum disordered state is expected. We show in Fig. the 
region corresponding to a reduced weight on the 32-site 
cluster. 

The reduction of Z is the first signal that the spinon 
moves away from the original site next to the vacancy 
at an average distance £ C onf to be determined. A quan- 
titative measure of this effect is provided by a careful 
inspection of the average local spin density (Sf ) around 
the vacancy in both the "bare" wavefunction and the 
true GS. Note that (Sf) in |^bare) gives the initial spin- 
spin correlation (SqS*) in the host GS. We start this 
analysis by examining the two extreme behaviors pro- 
vided by the Heisenberg model on the checkerboard and 
the kagome lattices reported in Fig. EJa) and (b) respec- 
tively. Clearly, the results for the checkerboard lattice 
show very short-ranged and incommensurate spin-spin 
correlations. In addition the spinon remains almost en- 
tirely confined on the N.N. site of the vacancy. In con- 
trast, on the kagome lattice, the spin-1/2 delocalizes on 
the whole lattice, a clear signature of deconfinement. Re- 
sults for the J1-J2-J3 model in Fig.[^a,b) for parameters 
corresponding to the two typical A and B points of the 
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FIG. 3: (Color on-line) Static hole (vacancy) spectral weight 
vs AF exchange parameters, (a) vs J2/J1 for J3 = and 
vs J3/J1 (for J2 = 0) as indicated on plot, (b) Along three 
different lines (J2 + Js)/Ji = est in the 2D ( J2/J1, J3/J1) 
parameter space. A and B refer to the points in the phase 
diagram of Fig. El 



phase diagram of Fig. (chosen because of a reduced Z 
factor) reveal completely new behaviors. First, we ob- 
serve for both A and B very short magnetic correlation 
lengths characterized by a fast oscillating decay (with 
the AF wavevector) of the correlations (see below). Note 
that no sign of incommensurability is seen in the oscil- 
lations of Fig. OJa) unlike in the classical spiral phase. 
Interestingly, the behavior of (Sf) in the "relaxed" |^ gs ) 
state differs drastically from the bare state with a much 
slower decay with distance [2l|. As seen in Fig. [SJ accu- 
rate fits can be realized by assuming a simple exponential 
decay together with an oscillatory behavior at wavevector 
(7r,7r) |22|. The very short correlation length £af, below 
one lattice spacing, is to be contrasted with the strikingly 
large confinement leng th £ C onf typically ranging from 2 to 
6 lattice spacings 23]. 

Let us now discuss some of the important implica- 
tions of such findings. First, we note that such non- 
trivial extended spin structure could be seen experimen- 
tally. Indeed, the substitution of a S=l/2 atom by a 
non-magnetic one (e.g. Zn 2+ for Cu 2+ ) which acts as 
a vacant site can be exactly described by our previous 
model. Moreover, the local spin densities (Sf) on the 
magnetic sites around a vacancy (spinless atom) in the 
bulk can be directly accessed by Nuclear Magnetic Reso- 
nance (NMR). It is important to notice that NMR would 
probe the spinon " wavefunction" in the "relaxed" state 
and not the host spin correlations. Note that newly de- 
veloped spin-polarized Scanning Tunnelling Microscopy 
(SP-STM) techniques might also allow to probe such 
atomic-scale spin structure around a vacancy on 
a surface. Secondly, at the theoretical level, the nu- 
merical evidence for a clear hierarchy of length scales, 
£conf ^> £af provides a strong argument in favor of the 
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FIG. 4: (Color on-line) Spin polarization in the vicinity of the 
vacancy (summed up on equivalent sites) for both the "bare" 
vacancy state and the GS for (a) the checkerboard lattice (32 
site cluster) and (b) the kagome lattice (30 site cluster). 
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FIG. 5: (Color on-line) Same as Fig. |5] (but for the modulus) 
for the J1-J2-J3 model with J2/J1 = 0.3 and J3/J1 = 0.2 
(a) and J2/J1 = 0.1 and J3/J1 = 0.4 (b) corresponding to 
points A and B in the phase diagram of Fig. El Fits using 
exponential forms are shown in dashed lines. The areas of the 
dots are proportional to the number of equivalent sites from 
the vacancy (entering in the fits). Dark and light symbols 
correspond to positive and negative values respectively. 



new class of DCP [l2j. Incidently, as seen in Fig. it 
is quite plausible that A and B lie indeed in the vicinity 
of the (supposed) Neel to VBS phase transition line, the 
paradigm of the DCP. 

Lastly, we examine the case of a mobile hole. This 
mimics an ARPES experiment in a Mott insulator where 
a single photo-induced hole is created or the case of a 
small chemical doping. The hole motion described as in 
a t-J model is characterized by a hole hopping amplitude 
t. For the unfrustrated t-J model, the hole dynamics 
has been successfully analyzed in term of holon-spinon 
boundstate [2^. Note that adding motion to the hole 
charge leads alone to a large reduction of Z, e.g. from 
0.93 (t = 0) down to 0.36 (Ji/t = 0.4) at small frustra- 
tion (J3/J1 = 0.05). As seen in Fig. as frustration 
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FIG. 6: (Color on-line) Four lowest energy poles of the sin- 
gle hole Green function (for hole momentum k — (jr, 0)) vs 
increasing frustration J3/J1. The GS energy of the undoped 
AF sets the energy reference. We assume here J2 = and a 
fixed ratio of Ji/t = 0.4 (J 3 /t = 0.4) for J 3 < Ji (J 3 > Ji). 
The areas of dots are proportional to the spectral weights 
l(^ihLl^bare)| 2 , n = 1, 4. To set the scale, Z ~ 0.04 
(n = 1) for J 3 = 0.6. 

is increased, (i) the low-energy spectral weight decreases 
further and (ii) the quasiparticle peak (at the bottom of 
the spectrum) is rapidly redistributed on several poles. 
This remarkable behavior indicates a severe weakening 
of the binding between the two constituents or, equiv- 
alently, a rapid increase of the size of the holon-spinon 
boundstate. A spectral weight Z below 0.01 is typical 
in this maximally frustrated region (e.g. Z ~ 0.008 for 
J3/J1 — 0.5). The dynamic hole problem bein g in fact 
smoothly connected to the case of a static hole 26], our 
data at finite t bring then additional strong evidence for 
the proximity of a DCP. 

To conclude, the confinement of a spinon liberated by 
introducing a vacant site or a mobile hole has been stud- 
ied in various frustrated Heisenberg AF. In the region of 
large frustration of the J1 — J2 — J3 model, an intermediate 
behavior between a strong confinement (as in the checker- 
board Heisenberg model) and a complete deconfinement 
(as on the kagome lattice) is observed, suggesting the 
emergence of a new length scale related to the confine- 
ment of the spinon. Its large value compared to the spin- 
spin correlation length supports the field-theoretic "de- 
confined critical point" scenario |l2| for the Neel-VBS 
transition. Furthermore, an interesting connection be- 
tween this real-space picture and features in the hole 
spectral function is established. 
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